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Abstract 

We consider here the vacuum structure in QCD with both quark and gluon 
condensates and a variational ansatz for the ground state. The method is 
nonperturbative using only equal time algebra for the field operators. We 
then find that a constrained energy minimisation of the Hamiltonian leads to 
a QCD vacuum with both quark and gluon condensates for a s > a c = 0.62. 
Pion decay constant and the charge radius of the pion seem to fix the QCD 
coupling constant a s as 1.28. The approach opens up possibilities of relating 
the mysterious vacuum structure with common place hadronic properties. 



1 



Typeset using REVTf^X 



I. INTRODUCTION 



It is now believed that quantum chromodynamics (QCD) is the correct theory for strong 
interaction physics of quarks and gluons. At low energies, however, the coupling constant in 
QCD becomes large leading to breakdown of the perturbative calculations. In this regime, 
the vacuum structure is also known to be nontrivial M with nonzero expectation values for 
quark and gluon condensates [EJ. Instability of QCD vacuum with constant chromomagnetic 
field or with vortex condensate formation has been studied since quite some time with 
a semiclassical approach M. QCD vacuum has also been studied with gluon or glueball 
condensates [§]|| as well as with nonperturbative solutions to Schwinger Dyson equations 
||. Further, a nontrivial vacuum structure with quark condensates in Nambu Jona Lasinio 
type of models @ has been seen to be consistent with low energy hadron physics. It is 
therefore desirable to examine the vacuum stucture in QCD with both quark and gluon 
condensates. 

We had proposed earlier a variational method which is nonperturbative with an explicit 
structure for the QCD vacuum. This has been applied to the case of gluon condensates for 
vacuum structure in SU(3) Yang-Mills fields demonstrating the instability of perturbative 
vacuum when coupling is greater than a critical value |J. Same methods have been applied 
to study the vacuum structure with quark condensates in QCD motivated phenomenological 
effective potential models |§ and in Nambu- Jona-Lasinio model [10|. For the ground state 
or vacuum this has been achieved through a minimisation of energy density, free energy 
density or thermodynamic potential depending on the physical situation. For the chiral 
symmetry breaking we had also considered a simple ansatz of taking the perturbative quarks 
having a phenomenological gaussian distribution in the nonperturbative vacuum [[II] . This 
appeared to describe a host of low energy hadronic properties as being related to the vacuum 
structure associated with chiral symmetry breaking. However, here no energy minimisation 
has been attempted. In the present paper we shall analyse the vacuum structure in QCD 
with both quark and gluon condensates and discuss its stability as opposed to taking only 



2 



gluon condensates [@-fj|| or only quark condensates [[7|P,|10| . 

As earlier, we shall take specific forms of condensate functions for quarks and gluons 
to describe a trial state for destabilised vacuum. Such an ansatz necessarily has limited 
dynamics. We shall circumvent this partially with the constraints that the SVZ parameter 
IT < G a ^ u G ailu > and the pion decay constant c n shall be correctly reproduced as experi- 
mentally known. The condensate functions still contain two parameters, over which energy 
is minimised. We then note that for a s greater than a critical coupling a c , perturbative 
vacuum destabilises with nonvanishing condensates in both quark and gluon sectors. 

We organise the paper as follows. In section II, we briefly recapitulate quantisation 
of QCD in Coulomb gauge and give an explicit construct for the nonperbutative vacuum 
with quark and gluon condensates. The ansatz for the QCD vacuum is similar to BCS 
ansatz of Cooper pairs in the context of superconductivity. In section III we consider the 
stability of such a trial state through an energy minimisation where pion decay constant and 
SVZ parameter are taken as constraints as above and discuss the results. In section IV we 
summarise our conclusions. We also notice that for a s ~ 1.28, the ansatz function of Ref. 
[ n | along with the correct charge radius of the pion is reproduced. 

The method considered here is nonperturbative as we shall be using the equal time 
quantum algebra for the field operators, but is limited by the choice of the ansatz functions. 
The technique has been applied earlier to solvable cases |L2| to examine its reliability and 



to ground state structure for electroweak symmetry breaking and cosmic rays |H| and to 



nuclear matter, deuteron or quark stars Jl4 



II. VACUUM IN QCD WITH QUARK AND GLUON CONDENSATES 

The QCD Lagrangian is given as 

(-'gauge (-"matter (-'inti (1) 

where 
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c 



^u 9 e = --G a ^{d^W\ - d u W% + gf ahc W\w:) + -G\ u G a ^, 



(2a) 



^matter = ^7^)^ 



(2b) 



and 



(2c) 



where W a ^ are the SU(3) colour gauge fields. We shall quantise in Coulomb gauge |1| and 
write the electric field G a oi in terms of the transverse and longitudinal parts as 



G a Qi + dif a , 



(3) 



where f a is to be determined. We take at time t=0 [§ 

W\{x) = (2vr)- 3 / 2 



dk -> t 

'a a i(k) +a a i(-k) )exp(ik.x) 



2u{k) 



and 



T G a 0l (x) = (2vr)- 3/2 2 / dk 



\| ^Y-(-a a i{k) + a a t (-kf) exp(ifc.f), 



where, u;(A;) is arbitrary ]T3| and for equal time algebra we have 



a a i {k),a b j {k'f] = S^Ayfysfi-k'), 



with 



The equal time quantization condition for the fermionic sector is given as 



(4a) 



(4b) 



(5) 



(6) 



(7) 



where i and j refer to the colour and flavour indices 
expansion for fermion field i/j at time t=0 given as [RBI 



We now also have the field 
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(2tt) 3 / 2 

where U and V are given by [ffC 



U r (k)c\ r (k) + V s (-k)5 l Is (-k) e ik - x dk, 



(8) 



- l( l\ 

U r (k) = -= 



U Jr \ V s (-k) = -j= 



<a.k\ 



V 1 / 



Vis 



(9) 



for free chiral fields. The above are consistent with the equal time anticommutation condi- 



tions with [|16 



4(^,4(^)1+ = 5 rs 5^8(k-k') = [eU$,4(#)t] 



(10) 



In Coulomb gauge, the expression for the Hamiltonian density, T 00 from equation (1) is 
given as 
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T 00 = : - T G a 0l T G a 0l + -W a i(-V 2 )W a i + gf ahc WtW)d l W ( r 

,2 



+ 9 -f abc r f w^w^ + I^rxw 



(11) 



where : : denotes the normal ordering with respect to the perturbative vacuum, say | vac >, 
defined through a a i{k) \ vac >= 0, c l Ir {k) \ vac >= and c l Ir {k)^ \ vac >= 0. In order to 
solve for the operator f a , we first note that 



/« = -VT - 9 f abc (V 2 ) <W C ). 



(12) 



Proceeding as earlier M with a mean field type of approximation we obtain, 

y 2 W\(x)+g 2 r hc f cde < vac' | ^(^(vT 1 ^-^) | vac' > d 3 W e (x)) 

= J a o(x), (13) 



where, 



r Q = g r bc wfo c oi - #7° 



(14) 
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As noted earlier [§,0], the solution of equation (13) does not suffer from Gribov ambiguity 
Clearly the solution for Wq(x) will depend on the ansatz for the ground state \vac' >. 
We shall now consider a trial state with gluon as well as quark condensates. We thus 



explicitly take the ansatz for the above state as H-ll 



\vac' >= U G Up\vac >, 



(15) 



obtained through the unitary operators U G and Uf on the perturbative vacuum. For the 
gluon sector, we have [§] 



U G = exp (B G i-B G ), 
with the gluon condensate creation operator B G ^ as given by 



(16) 



Bg ] = ^ / f{k)a%kya\{-kydk, 



(17) 



where f{k) describes the vacuum structure with gluon condensates. For fermionic sector we 
have, 



U F = exp (5 F f - B F ) 



(18) 



with HU 



BJ 



h(ky i(k) \ (a ■ k)c\{-k) 



dk. 



(19) 



Here h(k) is a trial function associated with quark antiquark condensates. We shall minimise 
the energy density for \vac > to analyse vacuum stability. For this purpose we first note 
that with the above transformation the operators corresponding to | vac' > are related to 
the operators corresponding to \vac > through the Bogoliubov transformations 



bt(k) 
bf(-kV 
for the gluon sector and 



coshf(k) —sinhf(k) 
—sinhf(k) coshf(k) 



\ ( am ^ 



<(-/c)t 



(20) 
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/ 



cos(h(k)) —(a ■ k)sin(h(k)) 
(a ■ k)sin(h(k)) cos(h(k)) 



^ ( Cj(k) ^ 



^Cj(-fe) j 



(21) 



dj(Jfe) 
dj(-jfe) 

for the quark sector. 

Our job now is to evaluate the expectation value of T 00 with respect to | vac >. For 
evaluating the same, the following formulae will be useful. 



< vac' I: W a Ax)W b Ay) :| vac' >= 5 ab x (2tt)- 3 / dke^^ A i? -(jfe), 



(22) 



< vac' | : T G a m {x) T GUy) : | vac' > = 5 ab x (2vr)- 3 / dM L{s '^ ^ i:j {k)u{k)F^{k). (23) 



In the above F±(k) are given as 



F±(k) = sinh 2 /(Jfc) ± 



sinh2/(fc) 



Similarly, for the quark fields we have the parallel equations given as 



where 



<: C(^)V^) :>w= (2tt)- 3 ^| (A^)^-*^^ 
<: i>i(x)4(^ :> vac ,= (2;;) V' /' (A^*))^*"^, 



A±(jfe) = ±-(7°sin2/i(fc) - 2(5 ■ fc) sin 2 /i(% 
2 



(24) 



(25a) 



(25b) 



(26) 



These relations will be used to evaluate the energy expectation value which is carried out 
in the next section. 

III. EXTREMIS ATION OF ENERGY FUNCTIONAL AND RESULTS 



We shall here proceed to evaluate the expectation value of the Hamiltonian of equa- 
tion (11). However for that we note that we have to know the auxilary field contribution 
(l/2)dif a dif a in equation (11) or equivalently the contribution arising out of the time like 



and the longitudinal components of the gauge field. As stated earlier we shall take a mean 
field type of approximation || to solve for W£ field as in equation (13). Then the solution 
for Wq field in equation (13) is given as 

J\$) 



W a (k) 



(27) 



k 2 + <j)(k) 

where Wq (k) and Jq (k) are Fourier transforms of Wq(x) and Jq(x) and, 4>(k) given through 



I/IN 3g 2 r dk' n ,,,./, , /2 (k 2 - k' 2 ) 2 , fc + A;' 



(28) 



Stt 2 J oo(k') --^^-'V" ' " 
Using equations QTI|), (p2|), ( |23| ) and (25), we then obtain the expectation value of T 00 
with respect to | vac > as 



f , __nn I / 

eo = < vac |: T :| t>ac > 
= CV + Ci + C*2 + C3 2 + C4, 



(29) 



where 



12iV / 
(2^)3 



dk\k\ sin 2 /i(/c), 



(30a) 



7T- 



oo{k)k 2 F_(k) dk, 



(30b) 



C 2 = <: l^M-V')^ :> 



A- 4 



7r 2 J u;(/c) 



F+(fc) dfc 



Cs 2 = <■■ ^f abe f aef W b i W e j W e i W f j :> vad 



'2g r k 2 
ix 2 J Lu(k) 



F + {k) dk 



(30c) 



(30d) 



and 



In the above, 



G x {k) = 3g 2 J dqF + (\q\) F_(|* + q\) 



"(I k + q\) 

"(\q\) 



x (i + (31a) 
and, the contribution from quarks, 

G a (jb) = -^a 2 J dq[sm2h( y q) sm2h(k + q) 

4 g • (fc + g) gin 2 w hk gin 2 + ^ (31b) 

Ml* +51 

and 0(fc) as given earlier in equation (28). As may be noted here, the contributions from 
the quark condensates to the energy density comes through auxiliary equation through Wq 
as well as from the quark kinetic term. 

We shall now minimise the energy functional eo of equation (29). For the same we shall 
take u(k) to be of the free field form with an effective mass parameter for the gluon fields 
given as 



cj(k) = ^k 2 + m 2 G . (32) 

Here the gluon mass parameter mg given through the self consistency condition || 

2g 2 



F + {k) dk (33) 



G vr 2 J u{k) 

arising from the sum of the single contractions of the quartic gluon field interaction terms 



of T 00 in equation (11) @|Tj 

The condensate functions f(k) and h(k) are to be determined such that the energy 
density eo in (29) is a minimum. In simple cases it could be possible to solve for the 



condensate functions through functional differentiation |10|JT2| . In the present case however 



the dependence of the energy density on the condensate functions is highly nonlinear and 
it is not possible to solve for the same through functional differentiation and equating it 
to zero. We therefore adopt here the alternative procedure of taking a reasonably simple 
ansatz for the condensate functions by parameterizing the same. We parameterize the gluon 
condensates as, with k — \k\, 

sinh/(fc) = Ae~ Bk2/ \ (34) 

which corresponds to taking a gaussian distribution for the perturbative gluons in the non- 
perturbative vacuum ||. Similarly for the function h(k) describing the quark antiquark 
condensates we take the ansatz 

A' 

tan2h(fc) = (eR2fc2 _ 1)1/2 - (35) 



The above is a generalisation of the ansatz of ref. [jTTJ which corresponds to A ' — 1 and 
vanishes when A ' — 0. It will be determined through energy minimisation. We note that 
the for free massive fermions tan2h(k) = m/\k\, so that the above ansatz corresponds to a 
momentum dependent mass given as 

m ( k ) = s„rJ A 1M/2 - ( 36 ) 



e 



1)1/2 



We may add here that such a definition of quark mass is the same as that obtained from 
the pole of the fermion propogator in a condensate vacuum |9||19[]. In the limit of zero 
momentum we then have the dynamically generated mass for the quarks given as 

A 1 

m q = — . (37) 
The relationship between the decay constant of pion and the quark condensate function has 



been discussed in ref. [11| and is given as 



A 7 f sin 2 2h(k)dk = C7r(m " )V2 (38) 



(2vr)3/2y v / 

where N n sin2/i(/c) is the wave function for the pion with 
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N~ 2 = / sin 2 2h(k)dk. 



(39) 



With the ansatz of equation (35) we then have 

/ V3 \ 2/3 



R 



7TC 7rA /m 



IT V " U 7T , 



(40) 



e x- + i-A' 

The above equation determines R as a function of A' when is known. This will be used 
when we extremise energy along with a parallel constraint for SVZ parameter as in equation 
(48). 

With the above ansatzes for the condensate functions the energy density eo may now be 
written in terms of the dimensionless quantities x = Rk and b = B/R 2 as 



where 



and 



with 



e (A) = ^f + h(A, b) + I 2 (A, b) + I 3 {A, bf + h(A, b)) 



3N f 



f / ™3 



x dx 



1 - 



.4 



'2 



1/2- 



e x -l + A' 2 

I x (A,b) = ^J uo(x,b)x 2 dx(A 2 e~ bx2 - Ae~ bx2/2 (l + A 2 e 
4 r x 4 dx 



2-bx 2 \l/2 



u(x, b) 



A 2 e- bx * +Ae- bx2/2 (l + A 2 e- bx ') 



2-bx 2 \l/2 



h(A, b) = 2 -{l ( A 2 e~ bx2 + Ae~ bx2 / 2 (1 + A 

it 2 J u[x, b) \ 



2 e -6x 2 U/2 



h(A,b) = 4 x (2tt)- 6 y dx 



G 1 (x) + G 2 (x) 
x 2 + d)(x) 



\2 e -bx'\l/2 



= 3/ J dx' [A 2 e- bx ' 2 + Ae~ bx ' 2 l 2 {\ + A 
x ^A 2 e- b ^ 2 + Ae- b ^ 2 / 2 (l + A 2 e- b ^ 2 ) l/2 ^j 



uj(\x + x'\) ( (x ' 2 + X ■ X ') 

X V ' ; - 17 X 1 + V 



u(x') 



X ' (x + X ') 



2 / ' 



(41) 

(42a) 
(42b) 
(42c) 
(42d) 

(42e) 



(43) 
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G 2 {x) = -^-g 2 j dx> 



A' 4 x' •(£ + £') 



{e x ' 2 - 1 + A' 2 )(e^+^ 2 - 1 + A' 2 ) \x '\\x + x '\ 



A' 2 \ 1/2 \/ ( A 12 \V2 

X 1- 1--^ r^r 1-1 



(e x ' 2 - I + A' 2 ) J J\ V ( e (^') 2 - 1 + A' 2 ) 



(44) 



and 



W= 3 6J&) { A2e ~ bX ' 2 + ^- te ' 2/2 (l + A 2 e- bx ' 2 f 2 



x(^-' 2 -^H^I)- («) 

In the above, 6) = (x 2 + fi' 2 ) 1 ^ 2 , with /i' = vhqR being the gluon mass. The self 
consistency condition M for gluon mass me in equation (33) yields 

I , 2 X ! d *j A 2 e-* 2 + Ae-* 2 / 2 (l + A 2 e-* 2 ) 1/2 ), (46) 



M 7T 2 J (X 2 + /i 2 )V2 v 

where \i = iiiq^/B. Clearly jj! and fi are related as fi' = fi/y/b. For given values of A and 
A', we first determine i? = R(A ') from equation (39), and then the parameter b from the 
SVZ parameter for gluon condensates through the equation 

£^ <: G%G a » u :> vac ,= 0.012 GeV\ (47) 

which explicitly gives || 

— j—^ x^(-h(A,b) + I 2 (A, b) + I 3 (A, b) 2 - h(A, b)) = 0.012 GeV\ (48) 

so that b is a function of A and A ' . 

We now minimise the energy density e with respect to the parameters A and A 1 with a 
self consistent determination for the gluon mass as in equation (47) for different values of the 
QCD coupling constant a s . We plot in curve 1 of Fig.l e as a function of a s . We note that for 
a s < 0.62 = a c the condensate functions vanish, so that for the present ansatz, perturbative 
vacuum is stable, whereas, for a s > a c , we have a transition to a nonperturbative vacuum. 
The energy density becomes more negative with increase of coupling and becomes about —55 
MeV/fm 3 when a s = 1.4. In curve 2 of the same figure we have plotted gluon condensate 
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parameter A min for different couplings. In contrast to ref. ||, near a s = a c this parameter 
has a discontinuity. In curve 3 we have plotted \J~B in units of fermi as a function of 
coupling. The length scale here is of the order of a fermi which appears to be reasonable 
in the context of confinement. In curve 4, we have plotted the gluon mass parameter mc 
in MeV as a function of coupling. This quantity appears to be constant and is around 300 
MeV as earlier [§], and is similar to the result of Cornwall || where we adjust for the QCD 



parameter A or, Monte Carlo simulations of lattice calculations |20|| . For such gluons QGP 
signals are considered for relativistic heavy ion collisions ||18|| . 

In Fig. 2 we have plotted different characteristics of the quark condensate functions. In 
curves 1 and 2 we have plotted the condensate parameters A' min and R respectively. They 
seem to increase as we approach a c from above. We have plotted in curve 3 the quark 
condensate value given as 



12 

(2V) 3 



< -ipip > 1/3 = — ^— / dksm2h(k) (49) 



1/3 



The condensation effect increases with coupling as expected. In curve 4 we have plotted the 
dynamically generated quark mass m q of equation ([37]). We may in fact identify the same 
as the parallel of constituent mass of quarks as obtained from chiral symmetry breaking. 

We may remark that for a s < a c , any nonzero trial functions make eo positive, which 
shows that for the present ansatz the perturbative vacuum is stable. 



IV. DISCUSSIONS 

We have considered here destabilisation of perurbative vacuum in QCD as a function of 
coupling constant through an explicit construction. We have taken an ansatz for nonpertur- 
bative vacuum with trial functions both for quark and gluon condensates. With the QCD 
Hamiltonian for quarks and gluons satisfying equal time algebra, we then do a minimisa- 
tion of the energy density with the constraints that the pion decay constant and the SVZ 
parameters are correctly reproduced. Minimisation of the energy density over remaining 
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parameters yields that vacuum with condensates is the preferred configuration when $ a s 
is greater than 0.62. In order to relate the vacuum structure of chiral symmetry breaking 
to c,,-, we exploit the results of Ref. [P-p"T] . 



It shall be amusing to consider the present results in the context of Ref. [JTTJ] where no 
extremisation was done. We find that for a s = 1.28, A' min ~ 1, and that the pion charge 



radius gets correctly reproduced. In fact, with A' = 1, we have [|11 



d2 p2 i p2 



(50) 



where 



N 2 

R\ = -E- I \^/sin2h{k)\ 2 dk 



(51) 



and 



R 2 



Nl 



sm 2 2h(k) k 2 R 4 tan 2 2h(k) 



4(1 -sin2h(k) 



dk 



(52) 



In the above is given in equation (^). With R ~ 0.96 fm, we then obtain that R c h — 0.65 
fm, which agrees with the experimental value of 0.66 fm [^TJ The new feature here is that 
we are able to relate the QCD coupling constant with low energy hadron properties. 

We may also note that the averge number of perturbative gluons or quarks in such a 
condensate vacuum is given as 



N G =< vac'la^fficiiffilvac 1 >= 2 x 8 x 



1 



2n ? ' 



sinh 2 f(k)dk. 



(53) 



and, 
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(2tt) e 



i? 3 7T 2 / ' J(h ' 



sin 2 h{k)dk 



1/2 



1 + A' 2 



1/2- 



(54) 



We then obtain that for example when a s = 1.28 corresponding to correct charge radius 
of pion, N G = 2A 2 j{jxBfl 2 ~ 0.233/fm 3 and that N q ~ 0.085/fm 3 . We also note that 
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for a s = 1.28, the "bag pressure" is given as (— cq) 1 / 4 = 140 MeV, which appears to be in 
general agreement with the phenomenological value of this parameter. 

We thus find that a constrained energy minimisation of QCD Hamiltonian can lead to 
both gluon and quark condensates with appropriate chiral symmetry breaking. Even we 
can have simultaneously both the charge radius and pion decay constant, features of chiral 
symmetry breaking, getting related to the QCD coupling constant. The method in fact 
seems to open up possibilities of linking up the (mysterious) vacuum structure of QCD to 
(common place) hadronic properties in an unexpected manner. 

We have focussed our attention on the structure problem. It will be desirable to include 
the dynamical effects in the light quark sector. Also, chiral symmetry breaking is important 
in the strange quark sector, and calculations adding condensates of the same can be carried 
out. Besides technical difficulty, this has the additional feature of Lagrangian mass being 
of the same order as dynamical mass and may need nontrivial extension of the technology. 
Work in this direction are in progress. 
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Figure Captions 

Fig.l: In curves 1,2,3 and 4, energy density e (in units of 100 MeV/fm 3 ), A min , \f~B (in 
units of fm) and effective gluon mass mc (in units of 200 MeV), are plotted respectively as 
functions of strong coupling constant a s . 

Fig. 2: In curves 1,2,3 and 4, A' min) R (in units of fm), the quark condensate (— < ipip >) 1 ^ 3 
(in units of 100 MeV) and dynamically generated quark mass m q (in units of 200 MeV), are 
plotted respectively as functions of the coupling constant. a s . 
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